June 28, 2011 0:13 



WSPC - Proceedings Trim Size: 9in x 6in grikin 



1 



o 

£^ ■ Kinetics of Chiral Phase Transitions in Quark Matter 

3 ■ 



(N 



Awaneesh Singh, Sanjay Puri 
School of Physical Sciences, Jawaharlal Nehru University, New Delhi - 110067, India. 

Hiranmaya Mishra 

O l' Theory Division, Physical Research Laboratory, Navrangpura, Ahmedabad - 380009, 

India. 
E-mail: hm@prl.res. in 



CD 



> 
o 

m 
in 



o 



- 1—1 

X 

5-H 



We study the kinetics of chiral transitions in quark matter using a micro- 
scopic framework (Nambu-Jona-Lasinio model) and a phcnomcnological model 
(Ginzburg-Landau free energy). We focus on the coarsening dynamics subse- 
quent to a quench from the massless quark phase to the massive quark phase. 
The morphology of the ordering system is characterized by the scaling of the 
order-parameter correlation function. The domain growth process obeys the 
Allen-Calm growth law, L(t) ~ t 1 / 2 . We also study the growth of bubbles of 
the stable massive phase from the metastable massless phase. 



The nature of the QCD phase diagram in the plane of temperature (T) 
and baryon chemical potential (/x) has been studied extensively over the 
last few decades. 1 ' 2 For fi = 0, finite-temperature perturbative QCD calcu- 
lations have been complemented by first-principle calculations like lattice 
QCD simulations. 3-5 However, for fi ^ 0, lattice QCD calculations are lim- 
ited to small values of /i. 6 In QCD with two massless quarks, the equilib- 
rium chiral phase transition is expected to be a second-order transition at 
zero baryon densities. However, in nature, the light quarks are not exactly 
massless and the second-order phase transition is replaced by an analytical 
crossover. This picture is consistent with lattice QCD simulations with a 
transition temperature T c ~ 140— 190 MeV. 7 Calculations based on differ- 
ent effective potentials, on the other hand, indicate the phase transition to 
be first-order at large [i and small T. This means that the phase diagram 
will have a tri-critical point, where the first-order chiral transition becomes 
second-order (for vanishing quark masses) or ends (for non-vanishing quark 
masses) at a critical end point in the phase diagram. The location of the 
tri-critical point (TCP)(/i £; , T E ) in the phase digram has been estimated 
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by lattice QCD simulations as fi E /T c ~ 0.94 and [i E /T E ~ 1.8. 8 

Heavy-ion collision experiments at high energies produce hot and dense 
strongly-interacting matter, and provide the opportunity to explore the 
phase diagram of QCD. While the high-T and small-/! region of the QCD 
phase diagram has been explored in recent experiments. Future heavy-ion 
collision experiments plan to explore the high baryon density regime, par- 
ticularly the region around the tricritical point. 9 The experiments at Rel- 
ativistic Heavy Ion Collider (RHIC) provide clear signals that the nuclear 
matter undergoes a phase transition to partonic phases at sufficiently large 
value of the energy density. However, the nature of the phase transition still 
remains an open question. We might mention here that lattice QCD ther- 
modynamics has an built in equilibrium assumption while the heavy ion 
experiments are essentially nonequilibrium processes. It is worthwhile to 
mention here that, in a phase transition process, information about which 
equilibrium phase has lowest free energy is not sufficient to discuss all possi- 
ble structures that the system can have. One has to understand the kinetics 
of the process by which the phase ordering or disordering proceeds and the 
nature of nonequilibrium structures that the system must go through on 
its way to reach equilibrium. 

It is not surprising, therefore that much attention has been focused 
on dynamics of chiral transition in dense quark matter, particularly near 
the critical point. The critical behavior and the fluctuations of conserved 
charges in the presence of spinodal decomposition in the context of chi- 
ral transition have been explored in Nambu-Jona-Lasinio (NJL) model. 10 
There have also been preliminary studies 11 of the kinetics of the chiral 
phase transition in a nonlinear sigma model coupled to quarks. In this con- 
text, the effect of dissipation in slowing down spinodal decomposition has 
also been studied in this model. 12 The dynamics of first-order phase tran- 
sitions has been considered recently in Ref. 13 Simulated by the findings of 
lattice QCD, few studies have also been considered for the case of smooth 
crossover, its dynamics as compared to a strongly first order transition. 14 
Further, a Langevin equation has been derived within NJL model describ- 
ing the behavior of the fluctuation above the critical temperature at finite 
chemical potential within a linear approximation. 15 

On the othcrhand, there has been intense research interest in the kinet- 
ics of phase transitions, and the phase ordering process that occurs after a 
rapid quench in system parameters e.g., temperature, pressure 16 ' 17 in vari- 
ous condensed matter systems. During the transition, the system develops 
a spatial structure of randomly-distributed domains which coarsen with 
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time. This domain growth process has been extensively studied in many 
condensed matter systems like ferromagnets, binary fluids, liquid crystals, 
etc. 

In the present work, we focus on the kinetics of the chiral phase tran- 
sition subsequent to a quench from the disordered phase (with vanishing 
quark condensate) to the ordered phase for quark matter. This means while 
the equilibrium configuration at zero temperature is an ordered state, the 
system finds itself in a configuration from the ensemble appropriate to a 
high temperature. We are thus interested in the far-from-equilibrium evo- 
lution of the system and pattern dynamics associated with it. 

To model chiral symmetry breaking in QCD, we use the two-flavor 
Nambu-Jona-Lasinio (NJL) model 18,19 with the Hamiltonian 

H = V> W ' V + 7 V) i> ta - G [{^f - (^7 5 tV0 2 ] • (1) 

i. a 

Here, m t is the current quark mass - we take this to be the same (ra; L = 
m) for both u and d quarks. The parameter G denotes the quark-quark 
interaction strength. Further, t is the Pauli matrix acting in flavor space. 
The quark operator ip has two indices i and a, denoting theflavor and 
color indices, respectively. This model exhibits a second-order chiral phase 
transition for massless quarks at small fi and high T, and a first-order 
transition at large [i and small T. 

To describe the ground state, we take an ansatz with quark-antiquark 
condensates: 20 

\vac) = cxp 

Here, , q are two-component quark and antiquark creation operators, and 
|0) is the perturbative chiral vacuum. Further, is a variational function 
related to the quark-antiquark condensate as 

2 

(^) = -A V [dksm[Mi(k)]. (3) 

This flavor-dependent function can be determined by minimizing the energy 
at T — 0, or the thermodynamic potential at nonzero T and density. In the 
mean-field approximation and, near the chiral phase transition, we write 



dk qf{k)\ff ■ k)hi(k)qf(-k) - h.c. 



|0>. 



(2) 
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down the expression for the thermodynamic potential as 

12 



(2tt) 3 /3 



dk {in [i + e -/3(VP+M^-^)" 

I + e -/3(Vfe 2 + M2 +AI )" | 



In 



12 



M 2 



M 2 

4G' 



(4) 



Here, we have taken vanishing current quark mass, and introduce M = 
—2gp s with p s = (ipip) being the scalar density and g = G[l + 1/(4AT C )]. 
The details of the mean-field approximation are reported elsewhere. 21 




Fig. 1. (a) Phase diagram of the Nambu-Jona-Lasinio (NJL) model in the (^t, T)-plane 
for zero current quark mass. A line of first-order transitions (I, green online) meets a 
line of second-order transitions (II, blue online) at the tricritical point (tcp). We have 
(MtcpjTtcp) ~ (282.58,78) MeV. The dot-dashed lines Si and S2 denote the spinodals 
or metastability limits for the first-order transitions. The open symbols denote 4 combi- 
nations of (fi, T), chosen to represent qualitatively different shapes of the NJL potential. 
The cross denotes the point at which we quench the system for b < 0. (b) Plot of 
Q (M, /3, fi) from Eq. (4) as a function of M. The (fj,, T)-values are marked in (a). The 
solid lines superposed on the potentials correspond to the GL potential in Eq. (5) with 
a from Eq. (6), and b, d being fit parameters (cf. Table 1). 



The phase diagram for the chiral transition in the (/^, T)-plane resulting 
from Eq. (4) is shown in Fig. 1(a). For the numerical calculation of the 
thermodynamic potential, we have taken here a three-momentum ultravio- 
let cutoff A = 653.30 MeV, and the four-fcrmion coupling G = 5.0163 x 10" 6 
McV~ 2 . 22 With these parameters, the vacuum mass of quarks is M ~ 312 
MeV. At T = 0, a first-order transition takes place at /i ~ 326.321 MeV. 
For = 0, a second-order transition takes place at T ~ 190 MeV. The 
first-order line (I) meets the second-order line (II) at the tricritical point 
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(/itcp,Tt C p) ~ (282.58,78) MeV. The first-order transition is characterized 
by the existence of metastable phases. The limit of metastability is denoted 
by the dashed lines Si, S2 in Fig. 1(a), referred to as spinodal lines. 

Close to the phase boundary, the potential in Eq. (4) may be expanded 
as a Ginzburg-Landau (GL) potential in the order parameter M: 



ft (M) = ft (0) + ^M 2 + h -M A + ^M 6 + 0(M 8 ), 



(5) 



correct upto logarithmic corrections. 10,19 In the following, we consider the 
expansion of potential ft (M) upto the M 6 -term. This will prove adequate 
to recover the phase diagram in Fig. 1(a), as we see shortly. The first two 
coefficients in Eq. (5) can be obtained by comparison with Eq. (4) as 



ft(0) 



1 

2G 



4^ / dkk 2 (in [l + e -«*-f ) 



dk k 



+ ln 



+ 



1 + e 



]} 



1 + el 3 ^-^ 1 + eP( k +^) 



(6) 



We treat the higher coefficients as phenomenological parameters, which are 
obtained by fitting ft (M) in Eq. (5) to the integral expression for ft in Eq. 
(4). There are two free parameters in the microscopic theory (/j and T), so 
we consider the M 6 -GL potential with parameters b and d. For stability, 
we require d > 0. 

In Fig. 1(b), we plot ft (M) - ft (0) from Eq. (4) as a function of M. We 
show 4 combinations of ([i, T) as marked in Fig. 1(a), chosen to represent 
qualitatively different shapes of the potential. The solid lines superposed 
on the data sets in Fig. 1 (b) correspond to the GL potential in Eq. (5) with 
a from Eq. (6), and b, d being fit parameters (see Table 1). 
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Fig. 2. Phase diagram in (b/(dA 2 ), a/(<2A 4 )) -space for the GL-free energy in Eq. (5). 
The typical forms of the GL potential in various regions are shown in the figure. The 
open symbols denote the (/i, T)-values marked in Fig. 1(a). The cross denotes the point 
where we quench the system for b < 0, shows all different possible quenches. 

The extrema of the potential in Eq. (5) are determined by the gap 
equation f'(M) = aM + bM 3 + dM 5 = 0. The solutions are M = 0, 
and M± = (— b ± \Jb 2 — 4ad)/(2d). For b > 0, the transition is second- 
order, analogous to an M 4 -potential - the stationary points are M = (for 
a > 0) or M = 0, ±M + (for a < 0). For a < 0, the preferred equilibrium 
state is the one with massive quarks. For b < 0, the solutions of the gap 
equation are as follows: (i) M = for a > b 2 /(Ad), (ii) M = 0, ±M+, 
±M_ for b 2 /(Ad) > a > 0, and (hi) M = 0, ±Af+ for a < 0. A first-order 
transition takes place at a c = 3b 2 / (16d) with the order parameter jumping 
discontinuously from M = to M = ±M + = ±(3|6|/4d) 1 / 2 . The phase 
diagram for the GL potential is shown in Fig. 2. The tricritical point is 
located at 6 tcp = 0, a tcp = [cf. Fig. 1(a)]. The 4 combinations of (//, T)- 
values marked in Fig. 1(a) are identified with the same symbols in Fig. 
2. 

Next, we study dynamical problems in the context of the above free 
energy. Consider the dynamical environment of a heavy-ion collision. As 
long as the evolution is slow compared to the typical re-equilibration time, 
the order parameter field will be in local equilibrium. We consider a system 
which is rendered thermodynamically unstable by a rapid quench from the 
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massless phase to the massive phase in Figs. 1 or 2. The unstable disordered 
state evolves via the emergence and growth of domains rich in the preferred 
phase. 16,17 The coarsening system is inhomogeneous, and we account for 
this by including a surface tension term in the GL free energy: 

n[M] = Jdr |M 2 + b -M A + ^M 6 + y (VM) * . (7) 

In Eq. (7), Q[M] is a functional of the spatially-dependent order parameter 
M(r), and K measures the energy cost of spatial inhomogeneities, i.e., 
surface tension. 

The evolution of the system is described by the time-dependent 
Ginzburg-Landau (TDGL) equation: 

!««0~r^ + .«.>. 

which models the over-damped relaxational dynamics of M(r, t) to the min- 
imum of il [M]. 23 Here, T is the inverse damping coefficient, and 8(r,t) is 
the noise term satisfying the fluctuation-dissipation relation: (8 (r, t)) = 
and (6{P,t')6(r 7 ',t")} = 2TT5(P -r 7 ')5 (f - t"). We use the natural scales 
of order parameter, space and time to introduce dimensionless variables: 
M = M M' (Mo = VWH); (£ = v^TR); t = rt' (r = (rial)- 1 ); 

6 = (rlap^T 1 / 2 /^! 1 / 2 ) 6'. Dropping the primes, we obtain the dimension- 
less TDGL equation: 

■77- M (r, t) = -sgn (a) M - sgn {b) M 3 - AM 5 + V 2 M + 9 (r, t) , (9) 
at 

where A = \a\d/b 2 > 0. For T = 10 MeV, as fi takes values (in MeV) 311, 
321.75, 328 and 335, the corresponding values of A are 0.067, 0.14, 0.206 
and 0.255 from Table 1. 

The following evolutions are performed in the dimensionless units of 
length and time. To get them in physical units, one has to multiply them 
with the appropriate dimensional parameters. Two unknown parameters 
are K, the strength of the surface tension and the inverse damping co- 
efficient T. These parameters are obtained from the GL coefficients (see 
Table 1). Let us note that the surface energy can be calculated to be 
a = \[K{o?I <1 /b) J dz(dM/dz) 2 . The surface energy (a) for quark matter 
is poorly known and varies from 10-100 MeV/fm 2 at small temperatures. 24 
To get an idea about the length scale we can take a = 50 MeV/fm 2 . For 
T = 10 MeV, we can estimate y/K/a =1.6 fm. This factor has to be multi- 
plied with the dimensionless length scale to get the physical spatial length 
scale in units of fermi. Similarly for the inverse damping coefficient T, we 
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take it to be of the order 2T/s, where s is a quantity of order l. 25 This 
leads to t = 2.6 t' fm/s. 

In the following, we consider the phase transition kinetics for two dif- 
ferent quench possibilities. The hrst case corresponds to high T and low 
baryon density (/x), where the quenching is done through the second-order 
line (II) in Fig. 1(a) or Fig. 2. The second case corresponds to low T and 
high baryon density (/z), where the phase conversion process can probe the 
metastable region of the phase diagram. This can be achieved by quenching 
through the first-order line (I) in Fig. 1(a) or Fig. 2. Note that both these 
scenarios can be studied using Eq. (9) by choosing a fixed value of A and 
appropriate values of a, b. In our simulation, we have used A = 0.14. This 
value of A corresponds to, e.g., (p, T)= (231.6 MeV, 85 MeV) or (321.75 
McV, 10 MeV) in Fig. 1(a). 

First, we focus our attention to the ordering dynamics in the context 
of the phase diagram of Fig. 2 for the case of b > 0. This corresponds to 
the second order transition case relevant for low chemical potential regime 
in the context of the phase diagram of Fig. 1(a). For b > 0, the chiral 
transition occurs when a < 0. We solve Eq. (9) numerically using an Euler- 
discretization scheme with an isotropic Laplacian. We have implemented it 
on a d=3 lattice of size iV 3 (N = 256), with periodic boundary conditions 
in all directions. The dimensionless mesh sizes are Ax = 1.0 and At = 0.1, 
which satisfy the numerical stability condition. We have further confirmed 
that the spatial mesh size is sufficiently small to resolve the interface re- 
gion. 26 

Figure 3 shows the evolution of a disordered initial condition for Eq. 
(9) with b > and a < 0. It corresponds to a temperature quench through 
the second-order line (II) in Figs. 1(a) or 2. The initial state consists of 
small-amplitude thermal fluctuations about the massless phase M = 0. The 
system rapidly evolves into domains of the massive phase with M ~ M + 
and M ~ — M+. The interfaces between these massive domains correspond 
to M — 0. The interface evolution is shown in the snapshots (frames on 
left) of Fig. 3. The frames on the right show the interface structure in a 
cross-section of the snapshots. 

The characteristic length scale L(t) of the domains grows with time. 
The growth process is analogous to coarsening in the TDGL equation with 
an M 4 -potential, where the coarsening is driven by kinks with the equi- 
librium profile M{z) = tanh(±2:/v / 2). 16 ' 17 The order-parameter correlation 
function C(r, t) shows dynamical scaling C{r, t) = f(r/L). The scaling func- 
tion f(x) = (2/7r) sin _1 (e _a: ) has been calculated by Ohta et al. (OJK) 27 
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t = 20 




Fig. 3. Interface evolution after temperature quench through second-order line (II) in 
Figs. 1(a) or 2. The d=3 snapshots on the left show the interfaces (M = 0) at t = 20, 
100 (in units of r). They were obtained by numerically solving Eq. (9) as described in 
the text with a < 0, b > 0, A = 0.14. The noise amplitude was e = 0.008. The frames on 
the right show a cross-section of the snapshots at z = N/2. 

in the context of an ordering ferromagnet. In Fig. 4(a), we demonstrate 
that C(r, t) shows dynamical scaling for the evolution in Fig. 3 shows the 
scaling property. Further, the domain scale obeys the Allen-Cahn (AC) 
growth law, L{t) ~ i 1 / 217 [see Fig. 4(b)]. Typically, the interface velocity 
v <~ dL/dt ~ where L~ x measures the local curvature of the inter- 

face. This yields the AC growth law. The same growth law has also been 
obtained via a closed time path formalism of relativistic finite-temperature 
field theory applied to the NJL model. 28 

Next, let us consider the case with b < in Fig. 2. In this first- 
order chiral transition occurs for a < a c — 36 2 /(16g?) (or A < A c = 3/16). 
For a < 0, the potential has a double- well structure and the ordering dy- 
namics is equivalent to M 4 -theory, i.e., the domain growth scenario is sim- 
ilar to Figs. 3 and 4. We focus on a quench from the disordered state (with 
M = 0) to < a < a c or < A < A c , corresponding to a quench between 
the first-order line (I) and Si in Figs. 1(a) or 2 - the corresponding points 
are denoted by crosses. The massless state (M = 0) is now a metastable 
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Fig. 4. (a) Scaling of order-parameter correlation function, C(r, i) vs. r/L, for A = 0.14 
at four different times. The different data sets collaps onto a master curve. The statistical 
data is obtained as an average over 10 independent runs on 4096 2 lattices. The length 
scale L(t) is defined as the distance over which C(r, t) decays from 1 —J- 1/2. The solid line 
denotes the OJK result. 27 (b) Time-dependence of domain size, L(t) vs. t, for A = 0.14. 
The coarsening process obeys the Allen-Cahn (AC) growth law, L(t) ~ t 1 / 2 . 



state of the M 6 -potential. The chiral transition proceeds via the nucleation 
and growth of droplets of the preferred phase (M = ±M+). The nucleation 
results from large fluctuations in the initial condition or thermal fluctua- 
tions during the evolution. 

In Fig. 5, we show the nucleation and growth process. At early times 
(t = 400) , the system is primarily in the M = phase with small droplets of 
the preferred phase. These droplets grow in time and coalesce into domains. 
The subsequent coarsening of these domains is analogous to that in Figs. 
3-4, through the interfacial kinks are slightly flatter in the M ~ region. 
In the late stages of growth, there is no memory of the nucleation dynamics 
which characterized growth during the early stages. 

For b < and < a < a c , we have also studied the growth of sin- 
gle droplets of the preferred phase (M = +M + ) in a background of the 
metastable phase (M = 0) for A = 0.14 [see Fig. 6(a)]. We start with an 
initial configuration of a bubble of radius i?o > R c sucn that M(r) — M + 
for r < i?o and M(r) = for r > R , where R c is the critical size of the 
droplet. This configuration is evolved as per Eq. (9) with the noise term. 
The noise term is taken sufficiently large to overcome the metastable bar- 
rier. The droplets have a unique growth velocity v(X), which depends on 
the degree of undercooling (A): v — > as A — > A~. We have obtained v(X) 
by undertaking a phase-plane analysis of the traveling-wave solutions of Eq. 
(9) with 8 = 0. The droplet interface corresponds to a saddle connection 
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t=400 




Fig. 5. Interface evolution after shallow temperature quench through first-order line (I) 
in Figs. 1(a) or 2. The d=2 snapshots on the left show the interfaces (M = 0) at t = 400, 
4000. They were obtained by solving Eq. (9) with b < 0, a c > a > and A = 0.14. The 
frames on the right show the variation of the order parameter along the diagonal. Notice 
that the metastable patches (M = 0) at t = 400 are absent at later times. 



between the fixed points +M + — > or — M + — > 0. The details of this analy- 
sis will be presented later. 21 In Fig. 6(b), we plot numerical results for v(X) 
vs. A along with our theoretical result. 

Before concluding, we briefly discuss the inertial counterpart of the 
overdamped TDGL equation (8). The TDGL equation contains first-order 
time-derivatives and is not Lorentz invariant. This is reasonable as we are 
considering finite-temperature field theory, where there is no Lorentz in- 
variance. However, in some situations, the inertial terms can also play an 
important role in the evolution dynamics. The inertial TDGL equation has 
the following dimensionless form: 

f) 2 M BM 

-jjp- + = - s g n («) M - sgn (b) M 3 - AM 5 + V 2 M + 6 (f, t) , (10) 

where 7 measures the relative strengths of the damping and inertial terms. 
In principle, Eq. (10) can be obtained from a microscopic field-theoretic 
description of the noncquilibrium dynamics of the scalar field at finite tern- 
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Fig. 6. (a) Growth of droplet of the preferred phase (M = -\-M+) in a background 
of the metastable phase (M = 0) for b < 0, a c > a > and A = 0.14. We show the 
boundary of the droplet at four different times. The inner circle corresponds to a droplet 
at time t = 20. (b) Droplet velocity v vs. A. The circles refer to numerical data, while 
the solid line corresponds to the result from a phase-plane analysis. 21 

peratures. 29 We have also studied ordering dynamics for the inertial case, 
and will present details of our results later. 21 Here, we mention the main 
results of our study. The ordering dynamics in the inertial case is analo- 
gous to that in the overdamped case, except that nucleation does not have 
a significant effect even during the early stages of evolution for quenches to 
b < and < a < a c . The droplets grow very rapidly and merge to form a 
bicontinuous domain structure characteristic of late-stage domain growth. 
The domain growth law is again the AC law, L(t) ~ t x l 2 . 

To summarize: we have studied the kinetics of chiral phase transitions 
in QCD. In terms of the quark degrees of freedom, the phase diagram is 
obtained using the NJL model. An equivalent coarse-grained description is 
obtained from a M 6 -Ginzburg-Landau (GL) free energy. The chiral kinetics 
is modeled via the TDGL equation, and we consider both the overdamped 
and inertial cases. We study the ordering dynamics resulting from a sudden 
temperature quench through the first-order (I) or second-order (II) tran- 
sition lines in Figs. 1(a) or 2. For quenches through II and deep quenches 
through I, the massless phase is spontaneously unstable and evolves to the 
massive phase via spinodal decomposition. For shallow quenches through I, 
the massless phase is metastable and the phase transition proceeds via the 
nucleation and growth of droplets of the massive phase. The merger of these 
droplets results in late-stage domain growth analogous to that for the un- 
stable case. In all cases, the asymptotic growth process exhibits dynamical 
scaling, and the growth law is L(t) ~ t 1 ^ 2 . Given the dynamical universality 
of the processes involved, our results are of much under applicability than 
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the simple NJL Hamiltonian considered here. We hope that our results will 
motivate fresh experimental interest in hot and dense quark matter, and 
that our predictions will be subjected to experimental verification. 
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